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ABSTRACT 

This paper consists of two related parts: In the first part we derive an expression of the moment of 
inertia (MOI) of a neutron star as a function of observables from a hypothetical r-rnode gravitational 
wave detection. For a given r-rnode detection we show how the value of the MOI of a neutron star 
constrains the equation of state (EOS) of the matter in the core of the neutron star. Subsequently, 
for each candidate EOS, we derive a possible value of the saturation amplitude, a, of the r-rnode 
oscillations on the neutron star. Additionally, we argue that a r-mode detection will provide clues 
about the cooling rate mechanism of the neutron star. The above physics that can be derived from 
a hypothetical r-mode detection constitute our motivation for the second part of the paper. In that 
part we present a detection strategy to efficiently search for r-modes in gravitational-wave data. 
R-mode signals were injected into simulated noise colored with the advanced LIGO (aLIGO) and 
Einstein Telescope (ET) sensitivity curves. The r-mode waveforms used are those predicted by early 
theories based on a polytropic equation of state (EOS) neutron star matter (Owen et al. 1998). In 
our best case scenario (a of order lCU 1 ), the maximum detection distance when using the aLIGO 
sensitivity curve is ~ 1 Mpc (supernova event rate of 3-4 per century) while the maximum detection 
distance when using the ET sensitivity curve is ~ 10Mpc (supernova event rate of 1-2 per year). 

Our results suggest that if aLIGO takes data in 2015, it may be possible to set constraints for the 
EOS of the neutron star remnant of the Messier 82 supernova (SN2014J) that occurred in January 
2014. Depending on the r-mode parameters and detection sensitivity we may be able to determine an 
upper bound (in the worst case) on the r-mode saturation amplitude or (in the best case) determine 
its value based on the assumed EOS. 

Subject headings: aLIGO, neutron stars, equation of state, r-modes, gravitomagnetic waves, gravita¬ 
tional waves, supernova, CFS instability 


1. INTRODUCTION 

A hypothetical detection of r-mode gravitational 
radiation from newborn neutron stars could have at 
least three major implications in our understanding 
of neutron stars: (i) explanation of the low rotational 
frequencies of the observed neutron stars when compared 
to their possible rotational frequencies at birth, (ii) set 
constraints on the equation of state of the matter in 
the core of the neutron star and (iii) set upper bounds 
on a and settle the debate about the magnitude of 
the saturation amplitude of the r-mode mass current 
oscillations on neutron stars. 

After some preliminaries in section 2, the paper splits 
into two parts: the first part consists of sections 3-7 and 
the second part consists of sections 8-11. In the first 
part we present the physics of a neutron star that can 
be derived from a hypothetical r-mode detection: (i) 
the MOI of the neutron star, (ii) the EOS of the matter 
in the neutron star nucleus and (iii) the saturation 
amplitude (a) of the r-mode mass-current oscillations. 
These results constitute our motivation for the second 
part, to perform a sensitivity study on r-modes: using 
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(time-shifted) eLIGO data, recolored with aLIGO and 
ET sensitivity curves, we examine the distances at 
which our decision making algorithms can be sensitive 
to r-mode signals from newborn neutron stars with a 
False Alarm Rate (FAR) of 0.1% and a False Dismissal 
Rate (FDR) of 50%. 

The first part of the paper is organized as follows: 
in section 3 we discuss the power dependence on the 
gravitational-wave frequency and the saturation ampli¬ 
tude, a, of the r-mode oscillations. This result is used in 
section 4 to argue that among all possible r-mode sources 
newborn neutron stars are the most promising sources 
of detectable r-mode gravitational waves. In section 5, 
we present the time frame after a supernova explosion 
in which we would expect to detect r-mode gravitational 
waves. We end the first part of the paper with sections 
6 and 7, where we present the motivation for a r-mode 
gravitational-wave search. In section 6 we derive a 
relation between r-mode gravitational waves and the 
MOI of the neutron star, while in section 7 we show how 
this relation is used to constrain the EOS of the neutron 
star matter and subsequently set constraints on the sat¬ 
uration amplitude, a, of neutron star r-mode oscillations. 

The second part of the paper is organized as follows: 
in sections 8 and 9 we discuss the choice of parameters 
used to construct the waveforms needed to design the 
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sensitivity study. In section 8 we determine the range 
of possible values for the initial spindown frequency, 0 o , 
of the neutron star, while in section 9 we argue about 
the choice of the range of values for the parameter a. 
These 2 parameters determine the waveforms used in 
the sensitivity study presented in sections 10 and 11. In 
section 10 we present our choice of waveforms and plot 
several of them demonstrating the frequency spin-down 
dependence on a and f a , while in section 11 we present 
the method of waveform injection and recovery as well 
as the results we obtained on the detection distances for 
each waveform that was used. The paper ends with a 
discussion on the results from our sensitivity study as 
well as suggestions for future work that will follow in a 
second paper. 

2. PRELIMINARIES 

In the late 1990’s, the r-mode toroidal pulsations 
of a neutron star became very promising for gen¬ 
erating strong gravitational-wave signals due to the 
Chandrasekhar-Friedman-Schutz (CFS) instability they 
exhibit (Friedman & Schutz 1978a,b). R-modes of any 
harmonic, frequency and amplitude are subject to this 
instability at any angular velocity of the star (Andersson 
1998; Friedman & Morsink 1998). Therefore, even 
the smallest toroidal perturbations in the velocity of 
the neutron star mass currents will keep increasing 
in amplitude. Since the energy source of the r-mode 
oscillations is the rotational energy of the star, these 
small perturbations can eventually reach energy values 
of the order of the rotational energy of the neutron star. 
That would imply that all neutron stars are unstable 
and would contradict our observations. Therefore, some 
r-mode damping mechanism must be active resulting in 
the r-mode oscillation reaching a saturation amplitude, 
thus preserving the stability of the neutron star. 

In considering the saturation amplitude, its normal¬ 
ization is such that values of order 1 carry energy of the 
same order of magnitude as the total rotational energy 
of the neutron star. Some authors have introduced 
damping mechanisms that can cause saturation at 
r-mode oscillation amplitudes of order 10 -4 — 10~ 2 di¬ 
mensionless units (Bondarescu et al. 2009), while others 
have introduced mechanisms that cause saturation at 
amplitudes equal to or larger than 10” 1 (Alford et al. 
2012). On the other hand, some studies have examined 
the hypothesis that, under certain conditions, r-mode 
oscillations may get suppressed when neutron star 
models with a solid crust are considered (Lindblom et al. 
2000). In this case, the r-modes would be completely 
suppressed when the temperature of the neutron star 
drops below 10 8 °K. However, we are interested in the 
very early stages in the life of a newborn neutron star, 
when its core temperature is still around 10 9 — 10 10 °K. 
Therefore, this solid-core suppression, if it exists, does 
not come into play until much later on in the evolution 
of the neutron star. The uncertainty exemplified in this 
situation calls for the design of a search for a r-mode 
gravitational radiation from newborn neutron stars to 
discover their saturation amplitudes. 


was the outcome of numerous earlier studies e.g. (Brink 
et al. 2004) that discovered a coupling between the 
r-mode and other inertial modes. Bondarescu et al.’s 
numerical study showed that for initial amplitudes of 
order 10 , non-linear couplings of the r-mode with 

other (neighboring in frequency) inertial modes can 
result in saturation of the r-mode amplitude at values of 
order up to 10 -2 (in their most optimistic case scenario). 

If the assumptions in the Bondarescu et al. model 
do not hold and the r-mode does not saturate at 
their suggested amplitudes, a different mechanism may 
saturate the r-mode oscillations at higher amplitudes. 
Such possibilities were considered by Mark Alford et 
al., who worked on models where the r-mode oscilla¬ 
tions can reach saturation amplitudes 2 to 3 orders of 
magnitude larger than those predicted by Bondarescu 
et al. They explored mechanisms that can saturate 
the r-mode oscillations at amplitudes of order 10 _1 - 1 
(Alford et al. 2012). They discovered that non-linear 
bulk viscosity will cause r-modes to saturate at large 
amplitudes and spin down the neutron star in much 
shorter time scales than Bondarescu et al. had suggested. 

3. ENERGETICS 

In this section we present the expression for the 
r-mode waveforms that was derived in (Owen et al. 
1998) specifically for neutron star matter obeying a 
polytropic EOS. We then derive the power dependence 
on the r-mode gravitational wave frequency, and the 
saturation amplitude a. This is used in section 4 where 
we argue that newborn neutron stars are the most 
promising sources for r-mode gravitational waves. 


3.1. Frequency evolution and energy of the r-mode 
oscillation 

Though very simplistic, the Owen et al. ’98 model is 
still a very good approximation for the early stages of the 
neutron star spin-down (Owen et al. 1998). More com¬ 
plicated numerical methods have shown that a r-mode 
saturation amplitude a = ICR 2 can result in a spin-down 
whose energy loss can be detected as gravitational radi¬ 
ation by aLIGO (Bondarescu et al. 2009). When this 
saturation amplitude is used in the ’98 model, we see 
that there is a good agreement in the angular velocity 
evolution of the neutron star up to several months after 
the start of the neutron star spin-down. The evolution 
of the angular velocity f1 of the neutron star in the Owen 
et al. ’98 model is described by 

dfl 20 <jQ 

dt t G r 1-crQ 

with a = a 2 , where a is the amplitude at which the r- 
mode oscillation saturates, Q is a dimensionless equation 
of state dependent parameter that takes the value Q = 
9.4 x 10 -2 for the polytropic neutron star model and tqr 
is the gravitational radiation (e-folding) time scale and 
is given by 

i i / o 2 y +i 

tgr t G r \nGp J 


The Bondarescu, Teukolsky and Wasserman ’09 result 


( 2 ) 
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where G is the gravitational constant and p is the average 
density of the neutron star (Owen et al. 1998). Combin¬ 
ing (1) and (2) for l = 2, tqr = —3.3 s and neglecting 
terms with powers of crQ higher than 1 we get 

dn 2A q n 7 

( nGp ) 3 

where A = as . Integrating (3) we get the angular 
velocity evolution of the neutron star to be 


fi(t) 



1 

6 


( 4 ) 


where fl 0 is the angular velocity at birth. An extension of 
this model that includes magnetic braking effects, is de¬ 
scribed in (Ho & Lai 2000) and (Staff et al. 2012). From 
(4), we can derive the time evolution of the neutron star 
rotational frequency (with f a being the initial frequency) 
to be 


fit) 


fo 



12(2tt) 6 A Q 
3.3 



( 5 ) 


Dividing numerator and denominator of the right hand 
side by f Q we get 


fit) 


which results in 


1 



i 

7.3xl0 3 AQ A 6 
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= 1.1 x 10 _ 2 O | a | 2 



( 7 ) 

( 8 ) 


Teukolsky (red) vs Owen (blue) model with alpha=0.01 for a spinning down ns 



Fig. 1.— R-mode waveform (red) predicted by the best case 
scenario of Bondarescu et al. numerical simulations where the r- 
mode oscillations saturate at an amplitude a of order 10~ 2 . For 
comparison this is plotted on the same axes with the waveform 
described by the Owen et al. model with a = 10” 2 . The two 
waveforms are in a very close agreement during the early stages 
(up until ~ 10 6 s ) of the frequency evolution. This is a good 
indication that even though the Owen et al. model is very simplistic 
and deviates from complicated numerical simulations for the total 
duration of the spin-down, it may be very accurate for the early 
stages of the frequency spin-down. In our sensitivity study we used 
2.5 X 10 3 s waveforms from the Owen et al model. 


due to the CFS instability. Therefore, the r-mode sat¬ 
uration amplitude was assumed to take values of order 1. 

This is not the case in later simulations (Bondaresku, 
Teukolsky, Wasserman, 2009) where the r-mode oscilla¬ 
tions exhibit non-linear couplings with daughter modes 
and the saturation amplitude (a) reaches maximum 
values of order 10~ 3 - 10 -2 . The models that predict 
these saturation amplitudes are discussed in section 
9.1. Data from the Bondarescu et al. model have been 
obtained and a comparison of the two evolutions (both 
cases with a = 10~ 2 ) is given in Fig.l. 


and the time needed for the neutron star rotational fre¬ 
quency to evolve from some initial frequency f a to a fre¬ 
quency value f < f 0 (using /® / 6 ) is given by 



1kHz A 6 
—) ‘ 


( 0 ) 


The Owen et al ’98 model depends on two parame¬ 
ters: the initial angular velocity = 2nf 0 discussed 
in section 2 and the r-mode oscillation saturation 
amplitude a. Using the upper bound of the angular 
velocity H ns ~ 7.2 x 10 3 rads _1 (discussed in section 
8.1) we conclude that we can take f a — 1.1 x 10 3 Hz. 
Since the r-mode gravitational wave frequency, f gw , is 
related to the neutron star rotational frequency, f ns , 
by fgw = 4/3 f ns we conclude that r-mode gravitational 
waves from newborn neutron stars would have initial 
frequencies bounded above by f gw ~ 1.5 x 10 3 Hz. 


Estimating the value of the parameter a is not 
as easy as estimating f Q . During the decade follow¬ 
ing the ’98 paper there was no known mechanism 
that could stop the r-mode amplitude from growing 


To estimate the energy stored in a neutron star r-mode 
and the power radiated as r-mode gravitational waves we 
need to know the values of f gw and a. From Owen et 
al. ’98 (Owen et al. 1998) the energy of the l = 2 r-mode 
(with M = 1.4 solar masses, R = 12 km) is given by 

E r = 0.82 x 10 ~ 2 \a\ 2 MR 2 n 2 = 1.3 x 10 38 / 2 s |a| 2 (10) 
and using f gw = 4/3 f ns we find 

E r = 0.73 X 10 38 / 2 ,Ja| 2 (11) 

For f ns = 1.5 x 10 3 Hz and a = 1 we find that the 
energy of the r-mode of a newborn neutron star is of 
order 10 51 ergs. 


3.2. Power dependence on the r-mode radiation 
frequency and satuarion amplitude a 

The gravitational-wave amplitude h a of a signal with 
frequency f gw = f at a distance d and power radiated E 
is given by (Bondarescu et al. 2009) 
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(in inks units) or 




( 12 ) 


hl~ 6-3 x 1(T 37 (13) 

The gravitational-wave amplitude at a distance d, the 
gravitational-wave frequency and the r-mode oscillation 
amplitude a are related by (Ho & Lai 2000) 


h 


O 


1.5 x 10" 23 


lMpc\ 
~d~ ) 


f y 

1kHz ) 


a 


(14) 


Substituting (14) in (13), we get a relation between the r- 
rnode amplitude, the frequency of gravitational radiation 
and the power radiated in r-mode gravitational waves 


E& 3.5 x 10 19 / 8 |a| 2 (15) 

Since the gravitational radiation power is propor¬ 
tional to the 8 th power of frequency we expect a very 
rapid decrease of the signal power. From (15), we 
see that by the time the frequency drops to 80% of 
the initial frequency, the power will drop to 2.8 % of 
the initial power. Therefore, a 0.8 x 10 3 Hz signal, 
will have radiation power ~ 36 times less than the 
power of a 1.0 x 10 3 Hz signal and therefore, a sig¬ 
nal to noise ratio (SNR) 36 times smaller than the 
SNR of a signal at 1.0 x 10 3 Hz. This means that we 
are more interested in the initial stages of the spin-down. 


The duration of the signal when frequency drops from 
1.0 x 10 3 Hz to 0.8 x 10 3 Hz ranges from about 3.5 x 10 2 s 
for a = 1 to about 3.5 x 10 6 for a = 10~ 2 (derived using 
(9)). Because of its duration, such a signal is classified 
as a long-transient gravitational-wave signal. For this 
reason gravitational-wave searches for the r-mode signals 
are suitable for cross-correlation type analyses (Thrane 
et al. 2011). 


by (1), it is assumed that non-linear fluid forces are 
coupling the r-mode to other modes of the star in such 
a way that gravitational radiation is able to extract 
energy from the overall rotational energy of the star, 
not just from the r-mode. So during this phase it would 
not be appropriate to use (11) as the total reservoir of 
energy that is being radiated into gravitational waves. 


4. R-MODE SOURCES 

In this section we use the power dependence on / and 
a as expressed in (15), to argue that the newborn neu¬ 
tron stars are the most promising sources of detectable 
r-mode gravitational waves. We also discuss the relevant 
types of electromagnetic triggers as well as their event 
rates. 


4.1. Comparing newborn neutron stars to other sources 

Apart from newborn neutron stars there are sev¬ 
eral other sources of r-mode gravitational radiation: 
starquakes on isolated neutron stars (Duncan 1998), 
pulsar postglitch relaxation (Rezania & Jahan-Miri 
2000; Clemens & Rosen 2004), accreeting low mass 
x-ray neutron stars (LMXB) in a binary system (Bon- 
darescu et al. 2007; Levin 1999), and oscillations of 
the remnant star (delayed collapse) during the merging 
phase (before collapsing to a black hole) of a compact 
binary coalescence (CBC) (Baumgarte et al. 1996; 
Shibata & Uryu 2000). From (15) we see that, due to 
their high angular velocities newborn neutron stars will 
emit the most powerful r-mode gravitational radiation 
out of all of the above sources. Even if we consider 
all the CBC as r-mode radiation sources (ignoring 
the low likelihood of a delayed collapse) the event 
rate is much lower than the supernova event rate 
(within the detection distances). The LMXB, the 
starquakes on neutron stars and the pulsar glitches 
may be more frequent events but the gravitational 
radiation emitted by these sources is much weaker than 
that (predicted to be) emitted by newborn neutron stars. 


From (15), we see that for an initial frequency of 
1.0 x 10 3 Hz the r-mode gravitational-wave power ranges 
from order ~ 10 46 ergss _1 for a = 10“ 2 , to order 
~ 10 5O ergss _1 for a = 1. These values are comparable 
to the power of the f-mode gravitational radiation, 
which is the most energetic of all neutron star oscillation 
modes. Using the energy stored in the neutron star 
f-mode as given in (de Araujo et al. 2005), with order 
ranging from 10 48 ergs to 10 5 ° ergs and also using the 
damping times of order 1.0 x 10 2 s, (Ferrari et al. 2003) 
the gravitational-wave power of these modes is of order 
10 46 ergss _1 to 10 48 ergss _1 , thus concluding that 
r-mode gravitational radiation may be (or better, it 
used to be thought of) as powerful as the neutron star’s 
most powerful gravitational radiation. 

Note that (15) cannot be taken simply by taking the 
time derivative of (11). This is because expression (11) 
is just the energy of the r-mode, and does note include 
the rotational energy of the neutron star itself. In 
the saturated phase of the r-mode evolution described 


4.2. Electromagnetic Counterparts 

The r-mode search from newborn neutron stars 
depends on electromagnetic triggers from supernova 
type-I and type-II explosions. From (14) we see that 
a r-mode detection will give an estimate for the ratio 
a/d. Therefore, to extract any information about the 
magnitude of a it is necessary to know the distance 
to the source. Distances to type-I supernova can be 
calculated using the standard candle method with an 
error between 5 — 10% (Branch & Tammann 1992). 
Distances to type-II supernovae can be calculated using 
the expanding photosphere method giving an error 
of 10 — 15 % (Kirshner & Kwan 1974; Schmidt et al. 
1994). Furthermore, our sensitivity study results (table 
1 in section 11.1) show that aLIGO can be sensitive to 
r-mode signals only from newborn neutron stars within 
our local group of galaxies. Since distances to galaxies 
in our local group are already known a supernova 
explosion within our local group would automatically 
give information about the distance to the hypothetical 
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r-mode gravitational radiation source. 


The amplitude, A, of the r-mode mass current oscilla¬ 
tions evolves like 


4.3. Event rates within the detection distances 

Considering both supernova types together we expect 
2 — 3 galactic events per century (Dragicevich et al. 
1999), 1 event every 2-3 years at a distance of 5Mpc 
and 1-2 events per year at a distance of lOMpc (Ando 
et al. 2005). However, the local supernova rates seem 
to be higher around Milky Way: In the years from 
2002-2005 4 events were detected (1-2 expected) within 
4Mpc and 9 events were detected (4-8 expected) within 
lOMpc. This information suggests that in the Milky 
Way’s neighborhood the supernova rates are higher than 
the predicted ones, meaning that the supernova rates 
around Milky Way’s local volume are higher than the 
rates at distant typical volumes. This is a motivation 
to improve our detection algorithms thus increasing the 
detection distances and cover as much of our local group 
of galaxies as possible. 

The latest supernova (SN2014J) occurred in January 
of 2014 in the galaxy Messier 82 (M82) in the nearby 
group of galaxies M81 and it is a type-I supernova. 
This galaxy is at a distance of 3.5 Mpc from the Earth. 
Estimates show a supernova rate in our Local group of 
galaxies (up to a distance of 1.5 Mpc from Earth) of 
3-6 per century (Mannucci et al. 2008; Li & White 2008). 


5. TIME FRAME FOR A R-MODE GRAVITATIONAL WAVE 
DETECTION 

In this section we discuss the physics of the neutron 
star matter via the use of the r-mode instability window 
to identify the time-frame for a hypothetical r-mode de¬ 
tection. An electromagnetic trigger will not necessarily 
coincide with the emission of r-mode gravitational waves. 
Depending on the physics of the r-modes (saturation 
amplitude and dissipation forces) and also the cooling 
mechanism of the neutron star (that depends on the 
EOS of the neutron star matter) the r-mode emission 
may start from several minutes to up to a year after an 
electromagnetic trigger. 


5.1. The r-mode instability survival window 
(temperature versus angular velocity) 

In a hot, newly born neutron star, the main dissipa¬ 
tion mechanism is due to viscous forces (Andersson & 
Kokkotas 2001; Kokkotas & Andersson 2001). For the 
simplest neutron star models, two kinds of viscosity are 
normally considered: bulk viscosity and shear viscosity. 
At high temperatures (above 10 10o K), bulk viscosity 
is the dominant dissipation mechanism while shear 
viscosity dominates for temperatures below 10 6 °K. 
Hence, there is a window between 10 6 — 10 10 °K where 
the r-mode instability may be active. 

The instability window is defined by the temperature 
interval in which the dominant dissipation mechanism 
of the r-mode oscillation is the emission of (r-mode) 
gravitational radiation. 


A 


exp 


it 



(16) 


where Q is the angular velocity of the neutron star and 
t is the (e-folding) time scale of the r-mode oscillations. 
The above expression can be expanded as 


A ~ exp (iClt) exp ( —- ) exp ( —- ) exp ( —- ) (17) 

V T g J \ TsJ V T b) 

where g stands for dissipation due to gravitational-wave 
emission, s stands for dissipation due to shear viscosity 
and b stands for dissipation due to bulk viscosity. Equa¬ 
tions (16) and (17) imply that the imaginary part, 1/r, 
of the frequency is given by 



(18) 


Assuming the r-mode mass perturbations, Sv, have a 
time dependence that is given by (16) then the pertur¬ 
bation energy has a time dependence like 


5v5v* ~ exp ^- J (19) 

Therefore, the time derivative of the perturbation energy 
is given by 


dE _ 2 E 

dt t 


( 20 ) 


and the r-mode instability condition can be expressed by 


1 <0 or — + — + — < 0 (21) 

T T g T s T b 

The conditions given by (21) imply that the energy of the 
r-mode increases with time. Since t s and r b are always 
positive then r g has to be negative and smaller than t s 
and Tf, such that 


1 

Tn 



1 

Tb 


( 22 ) 


for the r-mode instability to take place. In general 1/r 
is given in terms of the angular velocity of the neutron 
star H and its temperature T (Owen et al. 1998). Thus, 
the equation 1/r = 0 defines a curve of Tl C riticai(T). 
Above this curve, the instability condition is satisfied, 
therefore, this Q C riticai{T) curve defines the window over 
which r-mode instability will occur as seen in Fig.2. 


Bulk viscosity arises because of the pressure and 
density variations due to the r-mode oscillation driving 
the fluid away from beta equilibrium. Neutrinos carry 
away the r-mode oscillation energy lost during this 
process. The balance between neutrino cooling and 
viscous heating plays an important role in the spin-down 
evolution of a neutron star (Bondarescu et al. 2009). 
Also, the bulk viscosity can be strongly affected by the 
presence of hyperons in the neutron star core. In the 
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Fig. 2.— Instability window of the r-mode: The neutron stars 
are born with temperatures of the order 10 11 — 10 12 °K. At high 
temperatures (above 10 10 °K) bulk viscosity is the dominant dis¬ 
sipation mechanism. Shear viscosity dominates for temperatures 
below 10® °K. Hence there is a window between 10® — 10 10 °K 
where the r-mode instability is active and can have an astrophys- 
ically significant role (Kokkotas & Andersson 2001). The dotted 
curve defines the boundary ft cr itical(T) of the instability region, 
the dashed curves (mostly hidden underneath the solid curves) rep¬ 
resent the steady state, heating = cooling, and the solid lines show 
the numerical solution of the evolution equations for two fiducial 
initial rotational frequencies ft = 0.811 k and ft = 0.211 ^. As 
shown above, the higher the a value is the sooner the initial cool¬ 
ing phase will merge with the corresponding steady state curve. 
The above image is taken from (Alford & Schwenzer 2012). 

presence of hyperons the bulk viscosity coefficient is 
stronger and therefore, hyperon bulk viscosity becomes 
relevant at lower temperatures (Kokkotas & Andersson 
2001 ). 

Shear viscosity (Andersson & Kokkotas 2001; Kokko¬ 
tas & Andersson 2001) is used to describe in a 
macroscopic manner the scattering events associated 
with momentum transfer. Neutron-neutron scattering 
provides the most important contribution in a neutron 
star. However, in cases the core of the neutron star 
becomes superfluid, electron-electron scattering is the 
dominant effect that contributes to the shear viscosity. 
Therefore, a superfluid neutron star core can affect the 
spin-down evolution. 

5.2. Detection window after a supernova 

A neutron star is born with core temperatures of 
~ 10 12 °K. Right after the catastrophic collapse, the 
formation of a protoneutron star is defined as the 
collapsed star with neutrinos trapped inside its core. 
About 1.0 x 10 2 s after the collapse these neutrinos 
are released leading to a dramatic temperature drop 
to 10 10o K (Lattimer & Prakash 2007). The transition 
from a neutrino-opaque protoneutron star core to a 
neutrino-transparent core, resulting in the neutrino 
burst observed during a supernova, signals the death of 
the protoneutron star and the birth of a neutron star. 

After the newborn neutron star enters the instability 
window, it undergoes an initial cooling. The spin-down 


starts after the steady state (cooling=heating) curve 
is reached. Depending on the saturation amplitude 
and the EOS, this can be reached at temperatures 
anywhere between 10 9 — 10 10 (Alford & Schwenzer 
2012). That means the spin-down may start right after 
the neutrino burst (1.0 x 10 2 s after the collapse of the 
star) or up until the neutron star cools down to 10 9 °K. 
This time depends on the cooling rate of the neutron star. 

The modified Urea process (MUP) can decrease the 
temperature of the neutron star from T = 10 10 °K to 
10 9 °K in about a year and down to 10 s °K within a 
million years. However, if the direct Urea process (DUP) 
is doing the cooling, the core temperature will drop to 
10 9 °K within minutes and then drop down to 10 8 °K 
within days (Yakovlev & Pethick 2004). The EOS 
determines which cooling evolution, T(t), the neutron 
star will follow. For a given EOS, the transition rate 
from slow to fast cooling occurs in a given mass range. 
This range is very narrow because of the sharp neutron 
core mass density threshold for the direct Urea process 
to take place (Yakovlev et al. 1999). 

The growth timescale of the r-mode after it enters the 
instability window is of the order of ~ 40 s. Therefore, 
within 10 minutes the amplitude may grow up to order 
of 1 if not saturated before that. As saturation is reached 
the r-mode gravitational radiation will equilibrate at 
a particular heating=cooling curve (Fig.2). How long 
after the birth of the neutron star will this take place? 
From the above paragraphs we see that, depending 
on the equation of state and the r-mode saturation 
amplitude, the spin-down of a newborn neutron star 
may start anywhere from a few minutes up to ~years 
after the neutron star is born. Therefore we conclude 
that if aLIGO takes data in 2015, it may be possible 
to set some constraints for the EOS of the neutron star 
remnant of the Messier 82 supernova (SN2014J) that 
occurred in January 2014. 

6. MOTIVATION FOR THE R-MODE SENSITIVITY STUDY: 

MOI OF A NEUTRON STAR. FROM A R-MODE 
GRAVITATIONAL WAVE DETECTION 

In this section we derive a generalization of the 
r-mode waveform (6) so that the frequency evolution 
is expressed in terms of any EOS. We then derive an 
expression for the MOI of the neutron star as a function 
of the r-mode gravitational-wave observables, E, f and 
/. Subsequently we show how the value of the MOI 
sets constraints on the possible EOS of the neutron 
star matter. This result is used in section 7 where we 
derive an expression for the saturation amplitude, a , 
as a function of the EOS. Using this expression, a set 
of possible values of a can be derived, thus setting an 
upper bound for a. 

6.1. Estimating the moment of inertia of a neutron star 
from a hypothetical r-mode detection 

In section 2 we have assumed a polytropic EOS for 
the neutron star matter. The constant p in (8) as well 
as the waveforms expressed in (6) are dependent on this 
EOS. Therefore, a hypothetical r-mode detection would 
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set upper bounds on the r-mode saturation amplitude a 
based on the assumed EOS. Using the same hypothetical 
detection two different EOS may result in two different 
upper bounds for a. 

In what follows, we derive an EOS-dependent and 
a-independent expression (37) or (38) (which turns out 
to be the moment of inertia of the neutron star) that 
can be evaluated from a r-mode detection. Thus we con¬ 
clude that a hypothetical r-mode detection can constrain 
the EOS of the neutron star matter and can also set 
upper bounds on the values of a for each permitted EOS. 


and 


J 2 = 


1 

MR 4 



p(r)r 6 dr 


( 31 ) 


1 = 


8tt 

3 MR 2 



(32) 


Substituting (23), (24) and (30) in (29) (for m = 2) re¬ 
sults in the EOS-dependent expression for the time evo¬ 
lution of the angular velocity of the neutron star 


In this section we seek a generalization of the spin- 
down formula in (Owen et al. 1998), so that no assump¬ 
tions are made for the equation of state (EOS) of the 
neutron star matter. The gravitational radiation time 
scale is given by 


— = -A m MR 2m n' 2m+2 J m 


where 


32t rG / (to — l) m A 2 / to + 2 \ 2to+2 

c 2m+3 y (2m. + 1)!! J \?n + ly 


(23) 


(24) 


J m. — 


MR 2 " 


p(r) 


,2m+2 


dr 


The energy E m of the r-mode is given by 
1 


E m = -cUf YMR'Jr, 


(25) 


(26) 


The rate of change of the energy of the r-mode oscillation 
is given by 


E m =- (27) 

Tg W 

and hence from (23), (26) and (27) we obtain 

E m = A m a 2 M 2 R 2rn+2 Cl 2rn+4 J 2 m (28) 

During the saturated non-linear phase of the evolution, 
when the r-mode oscillation has a fixed amplitude and 
energy, (28) gives the power of the mode being radiated 
as r-mode gravitational radiation. During this phase the 
main contribution to the r-mode is the m — 2 harmonic 
and the angular velocity of the star evolves according to 
(Owen et al. 1998) 


2f l / a 2 Q A 2Cla 2 Q 

7~gw \1 J T~gw 


(29) 


where the last step assumes that q 2 <1 and Q is defined 
by 


Q = 


3 J 2 
21 


(30) 


? 2 

Cl = -3A 2 a 2 MR 4 ^-Cl 7 (33) 

where A 2 = 7.7 x 10“'°. Using the result relating the 
angular velocity, w = 2irf, of the r-mode gravitational 
radiation to the angular velocity, Cl, of the neutron star 
lo = 4/3 x Cl => f = 2/37r x Cl, we get the gravitational 
radiation frequency evolution given by 


/ = -2.5 x 10- 65 Mi? 4 a 2 ^y-/ 7 kg _1 m- 4 s 5 (34) 

leading to the EOS-dependent r-mode radiation wave¬ 
form 


where 


1 



(35) 


p = 1.5 x 10' 64 Mi? 4 a 2 ^-kg" 1 m- 4 s 5 (36) 

and f 0 is the initial frequency of the neutron star 
spin-down. 


Using (28) for m = 2, and (33) we get the result 


E 


8n 


—- = —— / p(r)r 4 dr 

CICl 9 J 0 


(37) 


or, in terms of the gravitational-wave frequency, /, 


E 

Jf 

A hypothetical r-mode detection will provide us with 
estimates of the model parameters f 0 and p. From 
these we can estimate the r-mode gravitational wave 
frequency, /, as well as its derivative, /. Furthermore, 
from equation (12) we see that given the distance, d, 
to the source of the r-mode gravitational radiation, 
the frequency, /, of the gravitational wave and the 
gravitational wave strain, h, then the power, E, of 
the r-mode gravitational radiation can be estimated. 
Therefore, we conclude that given a hypothetical r-mode 
gravitational wave is detected, the integral on the right 
hand side of (38) that gives the moment of inertia of the 
neutron star can be estimated. 


r R 

27r 3 / p(r)r 4 dr (38) 

Jo 


where 
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7. CONSTRAINING THE EOS OF NEUTRON STAR 
MATTER AND THE R-MODE SATURATION 
AMPLITUDE 

In this section we plot the moment of inertia versus 
radius and the moment of inertia versus mass rep¬ 
resentations of 17 EOS. From a hypothetical r-mode 
detection we can estimate the moment of inertia of the 
neutron star and hence estimate the possible radius and 
mass values for the allowed EOS. Subsequent estimates 
of either the radius or mass of the neutron star would 
further constrain the EOS. For each possible EOS we 
can then set upper bounds on the r-mode saturation 
amplitude, a. We also discuss the clues an r-mode 
gravitational wave detection will provide about the 
cooling mechanism of the neutron star. 


7.1. Numerical solution to the TOV equations 

Given a pair (p(P),P c ) of an EOS and a central pres¬ 
sure of a neutron star we can integrate the Tolrnan- 
Oppenheimer-Volkov (TOV) equations 

dP G(p(P) + £)(m+^) 
dr r (r- (39) 

and 

dm , 

—— = 47rr 2 p (40) 

dr 

from r = 0 ..r' and P = P c ..P[r'). At every integral step 
we evaluate the following integrals 


In(r') = / p(P)r n dr with n=2,4 and 6 (41) 

Jo 

which can be used to re-write (32) and (31) like 

i= w» u{R) a, “ 1 * = wiji UR) <42) 

respectively, as well as (38) like 

jj = -2 n 3 h(R) (43) 

where R is the value of r' for which 

P(r') = 0 (44) 

Using equation (41) with n = 2 and r' = R we can 
express the total mass of a neutron star of radius R by 

M{R) = 4ttI 2 (R) (45) 


Moment Of Inertia vs Radius representations of EOS 



Fig. 3.— Moment of inertia versus radius representation of the 
EOS. The horizontal lines correspond to a range of values of the 
moment of inertia (estimated from a hypothetical r-mode detec¬ 
tion), in this case I\ = 2.2 and I 2 = 2.3. The corresponding radius 
range is shown in Fig.5. 

star matter can be thought of as: (a) a curve in the I — R 
plane (Fig.3), (b) a curve in the I — M plane (Fig.4) and 
(c) a curve in the M — R plane (Fig.5); all curves being 
parametrized by P c , with each P c value corresponding 
to a unique point on each curve. Figures 3, 4 and 5 were 
plotted by integrating the TOV equations using the 
Runge-Kutta method for a set of 17 equations of state. 

A hypothetical r-mode gravitational-wave detection, 
with r-mode waveform parameters f 0 and p, will give 
an estimate of the moment of inertia, / (oc I 4 ), of the 
neutron star and hence intersect the corresponding I — R 
and I — M curves. The points of intersection represent 
the set of all possible mass and radius values allowed by 
the estimated value of /. The mass and radius values 
are shown in (Fig.3) and (Fig.4). Using these figures we 
see that a subsequent measurement of either the mass or 
the radius of the neutron star, would pick a single I — R 
or I — M curve and hence fix a choice of the EOS. The 
pairs of all possible ( M i 7 Ri ) points (i corresponding to 
each intersection point in (Fig.3) and (Fig.4)) are shown 
in (Fig.5) which is the standard M — R representation 
of the EOS. 


Upon constraining the EOS to a few options we can 
substitute the corresponding values of /, J 2 and R as 
given by (42) and (44) respectively in (36) to get 




= 3.0 x 10" 66 


a 2 (h(R)f 

R 2 C h{R )) 


(46) 


Solving (39), (40) and (41) (for the given p(P)) while 
varying P c , results in the functions R = R(P C ), h{R, Pc) 
and h(R,P c ) where I 2 and I 4 are proportional to the 
neutron star mass and moment of inertia respectively. 
Using these results we can express the total mass, M, 
and the moment of inertia, /, of the neutron star as 
M(P C ) and J(P C ). Expressing the macroscopic observ¬ 
ables, I, R and M in terms of the central pressure, P c , 
demonstrates that the equation of state of the neutron 


Rearranging (46) we can express a as a function of the 
value of p estimated from the hypothetical r-mode de¬ 
tection and the allowed EOS as follows 


a = 5.8 x 10 32 


Ry/ ph(R) 
h(R) 


(47) 


Therefore, we can set constraints (for each allowed EOS) 
on the r-mode saturation amplitude. 
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Moment Of Inertia vs Mass representations of EOS 



Fig. 4. — Moment of inertia versus mass representation of the 
EOS. The horizontal lines correspond to a range of values of the 
moment of inertia (estimated from a hypothetical r-mode detec¬ 
tion), in this case I\ = 2.2 and I 2 = 2.3. The corresponding mass 
range is shown in Fig.5. 


Mass vs Radius representations of EOS 



Fig. 5. — Mass versus radius representation of the EOS. Each 
value of the moment of inertia determined by a hypothetical r-mode 
detection corresponds to a point on one (or more) M-R curves. A 
range of moment of inertia values would correspond to a part of 
the M-R curves. Depending on what EOS curves are cut by the 
horizontal lines in Fig.3 and Fig.4 constraints may be put on what 
the possible EOS are for the neutron star matter. 


7.2. EOS, a and the cooling rate of the neutron star 

After the protoneutron star core becomes transparent 
to neutrinos, the cooling mechanism (and hence the 
time at which the neutron star hits the heating=cooling 
curve) depends on the EOS of neutron star matter. For 
non-superfluid neutron stars there are two distinct and 
well known cooling mechanisms: fast neutrino cooling 
via the direct Urea process (DUP) and slow neutrino 
cooling via the modified Urea process (MUP). For a 
given EOS there is a density threshold, p t h, beyond 
which DUP can occur. For the given EOS, there is 
a total mass threshold, M t h , required for the central 


density p c to be equal to pth- Neutron stars of total 
mass below M t h can only cool through MUP. 

The cooling evolution splits into two universal cate¬ 
gories: the (slow) cooling for neutron stars of mass less 
than the threshold mass, M t h , and the (fast) cooling for 
neutron stars of mass greater than the threshold mass, 
Mth■ The cooling curves in each category are universal 
i.e. they are independent of the mass or EOS. 

Assume the neutron star is born with an angular 
velocity 0 o that stays approximately constant (ignor¬ 
ing magnetic braking) until the r-mode gravitational 
radiation starts spinning down the neutron star. A 
hypothetical r-mode detection will provide us with the 
value of Q 0 as well as with the value of p as shown in 
(46). This value of p together with the chosen EOS can 
be used to determine the saturation amplitude a which 
in turn can be used to determine the heating=cooling 
curve as shown in Fig.2. The horizontal line 0 = in¬ 
tersects the heating=cooling curve at T = T r and hence 
T r , the temperature at which the neutron star starts 
spinning down due to r-mode gravitational radiation, 
can also be estimated. 

The time, t r , at which the neutron star temperature 
is equal to T r can be estimated by the initial detection 
time of the r-mode gravitational radiation. Furthermore, 
the time, t 0 , at which the neutron star is born is given 
by the time at which the neutrino burst occurs and 
may also be detected. Since the temperature of the 
neutron star at birth is generically approximately equal 
to 10 10 °K, t 0 gives the time at which the neutron star 
has this temperature. 

The time elapsed between t r and t a provides the 
time taken for the neutron star temperature to drop 
from ~ 10 10o K to T r . Therefore, the neutron star 
cooling rate can be estimated. This cooling rate may be 
consistent either with DUP (total mass has to be above 
Mth ) or MUP (total mass has to be below M t h)- If the 
cooling mechanism (implied by the cooling rate) is not 
consistent with the mass of the neutron star (implied by 
the chosen EOS) that means the chosen EOS is not the 
correct one. 

Constraining the EOS of a neutron star using cooling 
mechanism considerations as well as the value of the 
MOI of the neutron star (as predicted by a hypothetical 
r-mode detection) together with equations (46) and (47) 
are the main results of the first part of the paper. These 
results demonstrate the importance of a hypothetical 
r-mode detection as a source of information about the 
physics of the state of matter in the interior of the 
neutron star as well as the neutron star mass-current 
oscillations. 

This information carried by r-mode gravitational 
waves constitutes the motivation for the second part of 
the paper: to perform a sensitivity study on r-mode sig¬ 
nals from newborn neutron stars. This study starts with 
section 8, where we first determine the characteristics of 
a set of representative waveforms that were injected on 
background detector noise and then recovered using a 
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clustering algorithm. 


8. WAVEFORM PARAMETERS: RANGE OF POSSIBLE 
VALUES OF f 0 

In the sections that follow we are performing a 
sensitivity study to examine the distances at which the 
aLIGO and ET detectors (combined with the prepro¬ 
cessing and decision making algorithms) will be sensitive 
for such signals. To perform the sensitivity study we 
need to choose a representative set of waveforms whose 
choice depends on the range of values of the waveform 
parameters, f a and a. To cover the whole spectrum 
of possibilities we need to determine the highest theo¬ 
retically possible rotational velocities of neutron stars 
at birth and also determine the highest theoretically 
possible values of a. In this section we determine the 
possible range of values for the rotational velocity of 
neutron stars at birth. 


8.1. Theoretical upper hound for the neutron star 
angular velocity at birth 

The observed neutron star angular velocities are 
up to 25 times less than the theoretical upper bound 
of angular velocities that neutron stars may possess 
at birth. Conservation of angular momentum of the 
progenitor star results in much higher angular velocities 
at birth than those our observations reveal. Therefore, 
there must be a mechanism that causes the neutron star 
to spin down to the observed angular velocities. 

A progenitor mass between 8 and 25 solar masses will 
result in a supernova with a neutron star as a remnant. 
Numerical evolutions (Heger et al. 2003) have shown that 
such progenitors reach angular velocities very close to the 
critical ones (called Kepler angular velocities) at which 
the stars would fall apart due to centrifugal forces. Three 
different approaches (general relativistic with and with¬ 
out differential rotation as well as numerical evolution of 
a pre-supernova stellar model) point towards an upper 
bound of 7.2 x 10 3 rads -1 for the critical angular veloc¬ 
ity of such neutron stars. In the Newtonian approach, 
this limit is obtained using 


v escape 


2 GM 
R 


(48) 


where G is the gravitational constant, M is the mass of 
the neutron star and R is the radius of the neutron star. 
Equation (48) gives 


^ Kepler (N) = 1.6 ^/nGp (49) 

while general relativistic corrections give 


n 


_ 2 

Kepler (gr) — ^ 


\JkGp 


(50) 


Assuming a solid neutron star without differential rota¬ 
tion of uniform average density p ns = 5.5 x 10 1 ' kgm -3 
we find that newborn neutron stars can reach angular 
velocities of Vt ns = 7.2 x 10 3 rads -1 without falling 
apart due to centrifugal forces. 


In general relativistic models with differential rotation 
(Dimmelmeier et al. 2008) the iron cores of supernova 
progenitor stars may rotate at angular velocities in the 
range of 0.45 — 13rads -1 (f lKe P ier (am — 6.7rads -1 for 
an iron core with mass M = 1.4 solar masses and ra¬ 
dius R = 1.0 x 10 3 km). Using conservation of angular 
momentum we get 


n 


ns — 


( M lc 
\M ns 



(51) 


From (Dimmelmeier et al. 2008) we find that the iron 
cores of supernova progenitors have masses in the range 
of 1.2 — 1.8 solar masses and radii from 1.0 x 10 3 km 
to 2.5 x 10 3 km. Substituting, Tl ic = l.Orads -1 , 
Mi C = M ns = 1.4 solar masses, i?,; c = 1.0 x 10 3 km and 
R ns = 12 km, we find that the angular velocity of a 
new-born neutron star is about 6.9 x 10 3 rads -1 . This is 
consistent with the 7.0 x 10 3 rads -1 upper limit found 
by Hashimoto in 1994 (Hashimoto et al. 1994). 


Simulating the angular momentum evolution of 
massive stars (8 — 25 solar masses) in their spheri¬ 
cally symmetric models Heger et al. presented the first 
rotating presupernova stellar models that included a 
1-dim prescription for angular momentum transport and 
centrifugal effects (Heger et al. 2000). These authors 
estimated the neutron star angular velocity by assuming 
that the total angular momentum contained in the 
pre-collapse iron core is conserved during collapse and 
supernova phases and is deposited completely in a 
rigidly rotating neutron star of radius equal to ~ 12 km 
and moment of inertia given by 0.35 M. ns R 2 (Ott et al. 
2006; Lattimer & Prakash 2001). In this way, they 
estimated an initial neutron star angular velocity of 
6.3 x 10 3 rads -1 for their fiducial 20 solar masses (zero 
age main sequence) stellar model. 


8.2. Possible mechanisms responsible for the neutron 
star spin-down 

According to the observational data, the angular 
velocities of neutron stars range between ~ l.Orads -1 
to ~ 1.0 x 10 4 rads -1 . Typical (generic pulsar) neutron 
stars have angular velocities from ~ 20rads -1 to 

~ 3.0 x 10 2 rads -1 . The highest of these is up to a 
factor of 24 less than the theoretical upper bound of 
the initial angular velocity of the newborn neutron star. 
Several suggested processes that could lead to an early 
spin-down of a newborn (proto)neutron star include 
(Ott et al. 2006): 

(a) angular momentum redistribution by global hydro- 
dynamic angular instabilities 

(b) r-modes and gravitational radiation back-reaction 

(c) rotation-powered explosions 

(d) viscous angular momentum transport due to convec¬ 
tion 

(e) neutrino viscosity or dissipation of shear energy 
stored in differential rotation 

(f) magneto-centrifugal winds 

(g) early magnetic dipole radiation 

(h) late-time fail-back and 
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(i) anisotropic neutrino emission 

It is also possible that slowly spinning cores are 
the natural end-product of stellar evolution. An¬ 
gular momentum transport via magnetic processes 
yields angular velocities equal to ~ 6.0 x 10 -2 rads -1 
(Spruit & Phinney 1998) that is too low to explain 
the observational data. The authors had to rely on 
subsequent spin-up by off-center birth kicks to obtain 
angular velocities of ~ 3.0rads -1 to match the obser¬ 
vations. Heger et al. followed a similar approach for 
magnetic angular momentum transport during stellar 
evolution to obtain neutron star angular velocities of 
~ 9.0 x 10 2 rads -1 (Heger et al. 2003). Having in 
mind that stellar evolution theories with rotation are 
still improving and given the uncertainties in iron core 
angular velocities and angular momentum profiles we 
should be very careful in accepting any spin-down 
mechanism as the final theory. 

In this paper we assumed that the neutron star 
is born with initial angular velocities of up to 
fi ns — 7.2 x 10 3 rads 1 and the neutron star has 
a temperature range that allows the r-mode gravi¬ 
tational instability to survive. Since the rotational 
energy of the neutron star is the energy source of the 
r-mode oscillations, the emission of r-mode gravitational 
radiation is assumed to be responsible for the decreasing 
rotational energy of the neutron star. 


9. WAVEFORM PARAMETERS: RANGE OF POSSIBLE 
VALUES OF a 

In this section we determine (by doing a thorough 
examination of the literature review) the highest theo¬ 
retically permitted value of a. The value of a has been 
the subject of many debates since the earliest work on 
r-mode gravitational waves. Here we present several 
contradictory results that were published during the last 
15 years about the maximum possible value of a. The 
still ongoing debate on the subject made our choice of 
range of values of a a little tricky. 


9.1. Historic perspective of the r-mode saturation 
mechanism 

The hypothesis that newborn neutron stars may spin 
down due to r-mode gravitational radiation was based 
on three major discoveries: 

(a) The existence of a secular instability in rotating 
Jacobian ellipsoids describing an increase in the angular 
velocity due to gravitational radiation (Chandrasekhar 
1970).' 

(b) The proof that non-axisymmetric modes of the form 
e lm ^ ( (j> being the azimuthal angle) for all rotating stars 
are driven towards instability by gravitational radiation 
reaction (Friedman & Schutz 1978a,b). This instability 
is known as the Clrandrasekhar-Friedman-Schutz (CFS) 
instability. 

(c) The discovery of non-radial low frequency oscillation 
modes on white dwarfs similar to the Rossby waves 
(hence the name r-modes) in the Earth’s oceans and 


atmosphere (Papaloizou & Pringle 1978). 

Though some work on r-mode oscillations on neutron 
stars was done during the 1980’s (Saio 1982; Wagoner 
1984) it wasn’t untill 1998 that it was demonstrated 
that neutron star r-mode oscillations satisfy the re¬ 
quirements for the CFS instability (Lindblom et al. 
1998). Lindblom et al. showed that r-modes belong 
to the group of non-axisymmetric modes that can be 
driven unstable due to gravitational radiation. This 
theory predicted that (assuming the r-mode oscillation 
amplitude grows sufficiently large) r-mode gravitational 
radiation (primarily in the m = 2 harmonic) would 
carry away most of the angular momentum of a rapidly 
rotating newborn neutron star in the form of r-mode 
gravitational radiation. 

In the early years after the discovery of the r-mode 
instability, initial scenarios assumed that the r-mode 
amplitude would grow to order unity (Lindblom et al. 
2002) before some unknown process would saturate the 
mode (Porquet & Dubau 1999; Andersson & Kokkotas 
2001). However, there were speculations that some 
non-linear hydrodynamics of the star might limit the 
growth of the r-mode to very small values. This could 
be the result of r-mode oscillations leaking energy 
(due to non-linear couplings) into other inertial modes 
faster than gravitational radiation reaction force could 
restore it. Since 1998 several studies on the saturation 
mechanism showed a variety of results. 

Early work (Lindblom et al. 2000) showed that the 
hypothesis of r-mode gravitational radiation cannot be 
applied to cold (T < 10 8 °K) neutron stars whose crust 
is not in a fluid state. It was shown that in neutron 
stars colder than 10 8 °K the crust might be perfectly 
solid and that would completely suppress the r-mode 
instability. However, the question about the r-mode 
amplitude growth and saturation on hot (T > 10 8 °K) 
neutron stars remained open. 

In the absence of any saturation mechanism the 
r-mode amplitude would grow to values of order 1. This 
idea was supported by early relativistic simulations 
that showed that on a fixed neutron star geometry 
no r-mode amplitude saturation occurs even at large 
amplitudes (Stergioulas & Font 2001). These findings 
were further reinforced by the results of fully non-linear, 
3-dim numerical simulations based on Newtonian hy¬ 
drodynamics and gravitation. This study investigated 
the growth of the r-mode and found that the r-mode 
amplitude reached values of order 1 (Lindblom et al. 
2001, 2002). However, both Newtonian and relativistic 
3-dim hydrodynamical simulations were (and still are) 
severely limited by computational time and cannot 
probe the timescales on which r-modes may saturate. To 
bypass this limitation Lindblom et al. used artificially 
enhanced gravitational radiation force (increased by a 
factor of 4500) in order to evolve the r-mode amplitude 
and allow it to reach values of order 1. 

Later on the problem was attacked analytically with 
the use of a weakly non-linear perturbation theory to 
study the non-linear interactions of the r-modes with 
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other inertial modes (Schenk et al. 2002; Arras et al. 
2003; Brink et al. 2004). Calculations on a network of 
~ 5000 inodes with about 1.3 million couplings showed 
that the r-mode amplitude saturated at small values 
of order 10^ 4 . Furthermore, they showed that under 
some circumstances, it is possible to show that the 
r-mode amplitude evolution is dominated by just one 
three-mode coupling (r-mode plus 2 daughter modes). 

The triplet of r-mode plus the 2 daughter modes, 
together with assumptions on the neutron star cooling 
mechanism and hyperon bulk viscosity, was used in 
further work on the r-mode instability: (a) describing 
the rotational frequency evolution of accreeting neutron 
stars (Bondarescu et al. 2007) and (b) model the 
spin-down of a newborn neutron star by the emission 
of r-mode gravitational radiation (Bondarescu et al. 
2009). The work on (a) showed that accreeting neutron 
stars undergo a cyclic thermal runaway evolution during 
which the r-mode amplitudes remain saturated at values 
of order 10~ 5 while work on (b) showed that the r-mode 
amplitudes in newborn neutron stars saturate at values 
of order 10 -4 to 10~ 2 . 

Using a single triplet of modes, the lowest para¬ 
metric instability threshold is very sensitive to the 
internal neutron star physics. The authors considered 
a model with a neutrino cooling mechanism, via a 
combination of fast and slow processes, viscous heating 
due to hyperon bulk viscosity and spin-down due to 
gravitational radiation and magnetic dipole radiation. 
Two necessary assumptions (that are also limitations) 
of this single-triplet model and were made so that 
r-mode amplitude saturation is achieved are: (a) high 
hyperon bulk viscosity and (b) fast neutrino cooling 
mechanism. The fast neutrino cooling assumes certain 
restrictions on the equation of state of the neutron star 
matter that may not be true for every neutron star. 
For example fast neutrino cooling by direct Urca-like 
processes requires exotic phases of matter in the core 
of neutron stars (Yakovlev & Pethick 2004; Yakovlev 
et al. 1999). Furthermore, if real neutron stars have 
low hyperon bulk viscosity the energy dissipated by the 
single triplet of modes used in (Bondarescu et al. 2007, 
2009) may not be sufficient to stop the growth of the 
r-mode amplitude. 

If the Bondarescu et al. assumptions for the neutron 
star model (based on the single triplet of modes) are 
not satisfied, several triplets of modes may be required 
to saturate the r-mode amplitude. In that case, several 
parametric instability thresholds would be passed before 
saturating the r-mode amplitude thus making it not 
clear at what order the latter would saturate. If the 
r-mode amplitude passes its first parametric instability 
threshold value, other near-resonant inertial modes 
are expected to be excited (via energy transfer from 
the r-mode) signaling the importance of non-linear 
effects (Bondarescu et al. 2009). This is a fundamental 
property of inertial modes which does not depend on the 
details of the equations of state or on a slow-rotation 
approximation. The r-mode frequencies reside inside 
an extended region where there is a sufficiently large 
number of inertial modes with similar frequencies, so it 


is easy for resonance conditions to be satisfied. 

Even if several multiple triplets of modes are included 
in simulating the non-linear effects, depending on how 
reliable the mode-coupling estimates of the r-mode 
saturation amplitude is, it may not be certain that the 
energy transfer to inertial modes will stop the insta¬ 
bility at low r-mode amplitudes. If the mode-coupling 
mechanism does not saturate the growth of the r-mode 
early enough Bondarescu et al. pointed out that other 
mechanisms like suprathermal bulk viscosity might then 
become relevant. The current scientific consensus is 
that the mode-couplings saturating the amplitude of the 
r-mode is a fundamental result and will remain valid 
even with an improved understanding. Nevertheless, 
it is thought very unlikely that the r-mode amplitude 
could grow larger than 10 _1 . 

In the absence of experimental evidence several 
authors still entertain the idea of r-modes saturating at 
amplitudes of order equal to or greater than 10” 1 . It was 
shown that if the non-linear coupling mechanism does 
not come into play the r-mode amplitude would keep ris¬ 
ing exponentially and eventually reach the suprathermal 
regime (where the chemical equilibrium becomes greater 
than KbT). In this high r-mode amplitude regime, 
the damping due to (suprathermal) bulk viscosity 
increases dramatically (non-linearly) with increasing 
r-mode amplitude. This causes the viscous damping to 
overcome the gravitational instability and saturate the r- 
mode at amplitudes of order 10 _1 — 1 (Alford et al. 2012). 

Considering the results found by the above studies it 
seems that the most likely values the r-mode amplitude 
may saturate at is anywhere between 10 -4 to 10 _1 . In 
our sensitivity study we considered r-mode waveforms 
of three saturation amplitudes 10” 3 , 10 -2 and 10 _1 . 
Considering the worst case scenarios (that r-mode 
amplitudes saturate at values less than 10 -2 ) our best 
chances to settle this open question is to design a 
targeted search of r-mode gravitational waves from 
newborn neutron stars. This requires some waiting time 
and some preparation in order to be ready for the next 
local group supernova. 

10. GRAVITATIONAL WAVE DETECTORS AND 
SIMULATED WAVEFORMS 

In this section we discuss the sensitivity level of the 
aLIGO and ET detectors We also present several plots of 
the frequency evolution dependence on f a and a (accord¬ 
ing to the Owen et al. ’98 model) and how these wave¬ 
forms are used to construct the signals that were injected 
into the time-shifted data. 

10.1. Sensitivity of the gravitational-wave 
interferometers 

The past decade has seen the development of a world¬ 
wide network of gravitational-wave interferometers. The 
initial LIGO detectors achieved design strain sensitivity 
of ~ 2 x 10 -23 Hz -1 / 2 in the most sensitive frequency 
range of 100 — 130 Hz. The next generation of interfer¬ 
ometers such as aLIGO and advanced Virgo (aVirgo) are 
expected to begin taking data starting in early Fall 2015. 
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The upgraded aLIGO/aVirgo experiments are expected 
to achieve a factor of ten higher strain sensitivity than 
initial LIGO/Virgo. In our simulation study below we 
show that aLIGO can probe r-rnode gravitational-wave 
spin-down signals to newborn neutron stars within the 
Local Group of galaxies. 



Fig. 6. — The sensitivity curves of aLIGO and ET: aLIGO sim¬ 
ulated noise shows a maximum sensitivity of ~ 5 X 10 -24 Hz —1 ' 2 
at 110 — 130 Hz while ET simulated noise shows a maximum sen¬ 
sitivity of ~ 2 X 10 _2S Hz -1 / 2 at 110 — 130 Hz. 

10.2. Choice of parameter values and waveform 
simulation 

The Owen et al. ’98 model was used to produce the 
waveforms injected in the MC and eLIGO noise colored 
with the aLIGO sensitivity curve and also the eLIGO 
noise colored with the ET sensitivity curve. In this 
model there are two parameters: the initial frequency 
fo and the saturation amplitude a. We assume that the 
gravitational-wave frequency of the signal is between 
600 — 1600 Hz (as shown in table 1). We also assume that 
the signal will be powerful enough for detection for at 
least 2500 s. This signal duration is enough for the power 
to drop to 17% of its initial value for waveforms with 
highest f 0 and a values. However, the power drops only 
to 87 % of the initial value for waveforms with the lowest 
f 0 and a values. In the latter case, longer durations 
would be required to increase the SNR of the signal 
(when clustering pixels in the detection algorithms) and 
hence increase the detection distance. However, using 
2500 s duration and 600 — 1600 Hz frequency span results 
in ft-maps of 2.5 x 10 6 pixels. To run the background 
and sensitivity study algorithms and to create ft-maps of 
this size requires memory that approaches the memory 
limitations of conventional CPUs. 

We used three f a values (700, 1100, 1500 Hz) and three 
a values (10 3 , 10 —2 , 10 -1 ) resulting in nine waveforms. 
The dependence of the waveforms on f a and a are shown 
in Fig.7 and Fig.8 


Owen model with a = 0.1 and 700 < f < 1600 
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Fig. 7.— Frequency evolution dependence on the initial frequency 
f Q . The higher the f Q value the steeper the spin-down. 


Owen model with f =1500 and 10 2 < a < 1 

0 



Time (sec) 


Fig. 8.— Frequency evolution dependence on the saturation am¬ 
plitude a. The higher the a value the steeper the spin-down. 

A gravitational wave emitted by a neutron star is a 
linear combination of the two polarizations, h + and h x 
and can be expanded in the detector’s frame like (Owen 
2010 ) 

h(t) = F + (t, if)h + (t) + F x (t, ip)h x (t) (52) 

where the coefficients F + (t,if) and F x (t,if) are the de¬ 
tector antenna pattern functions for the two polariza¬ 
tions, t is the time in the detector frame and i/’ is the 
polarization angle. The waveforms for the two polariza¬ 
tions are given by 

h + (t) = h 0 ^ + ^° S ^ cos$(t) (53) 

and 

h x (t) = h 0 cos (6) sin $(t) (54) 

where h 0 is the corresponding gravitational-wave strain 
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amplitude as measured on Earth, 9 is the inclination an¬ 
gle and <1>(i) is the time dependent phase of the gravita¬ 
tional wave in the detector’s frame and is given by 

$(f) = <f>(t 0 ) + = 4>(t 0 ) + 2nf gw (t l )dt l (55) 

where to is the fiducial start time of the observation and 
f gw is the gravitational wave frequency in the detector’s 
frame. The angle 9 can be set equal to zero by assuming 
that the direction of propagation of the gravitational 
wave is perpendicular to the plane of the interferometer. 


Owen spindown waveform with f = 1500Hz and a=0.1 



Time (sec) 

Fig. 9. — This is an example of a h+(t) gravitational-wave strain 
that was injected in simulated aLIGO and ET noise. 


11. INJECTION RECOVERY AND DETECTABILITY 

This section starts with a discussion on matched 
filtering, the optimal detection algorithm. This al¬ 
gorithm requires accurate knowledge of the signal we 
are searching for. Though this is not the case for the 
r-modes (due to the wide range of possibilities for the 
parameter values) matched filtering gives theoretical 
upper bounds and a measure of evaluating the efficiency 
of the clustering algorithm used in our present sensitivity 
study. Our study is designed to assess the detectabil¬ 
ity of r-mode gravitational wave signals injected on 
(time-shifted) data colored with the sensitivity curves of 
aLIGO and ET gravitational wave detectors. Focusing 
on realistic excess cross-power searches, we show that 
second-generation detectors such as advanced LIGO 
(aLIGO) may be able to probe neutron star spin-down 
models out to astrophysically interesting distances of 
~ 1 Mpc. Third generation detectors such as Einstein 
Telescope (ET) may be able to probe such signals out 
to distances of ~ 10 Mpc. 


11.1. Matched filtering estimates 

For gravitational-wave signals of a known form, the 
theoretically optimal search strategy is matched filtering 


(Owen & Sathyaprakash 1999). However, the r-mode 
signals we consider here are not suitable for a search 
with a matched filtering detection algorithm due to the 
lack of accurate models that are available to generate 
r-mode waveforms. A search with matched filtering 
must include a template bank that spans the space of 
all possible signals and unknown parameter values. An 
incomplete template bank can result in faulty upper 
limits if the true signal falls outside the template space. 
In order to generate a complete template bank, we 
require firm knowledge of the details of the waveform’s 
phase evolution. In presenting the above model, we have 
not aspired to this degree of accuracy. Moreover, even if 
a complete and accurate model could be written down, 
there may be computational challenges associated with 
performing the search, especially for long signals with 
many parameters. 

Nonetheless, it is useful to compare the detection 
distances calculated using the excess-power technique 
to estimates for what can be achieved with matched 
filtering as this places an upper limit on the detection 
distance that can be achieved through improvements 
to the data-analysis scheme. As before, we assume 
optimal orientation of the source, optimal orientation of 
the detector network, a false-alarm probability <0.1% 
and a false dismissal probability <50%. We find that 
highly idealized matched filtering allows us to extend 
the detection distance up to factors of 10 - 20 as shown 
in table 1. This factor varies depending in part on 
the efficiency of the pattern recognition algorithm for 
different signal types. 


11.2. Spectrograms and the seedless clustering detection 

algorithm 

To estimate the distances at which we can see the 
signals predicted by the Owen et al. ’98 model using 
various combinations of the model parameters (/ D ,a), 
we simulate an excess cross-power search for long-lived 
gravitational waves associated with a well-localized 
electromagnetic counterpart, in the r-mode case this 
would be a supernova explosion. This simplifies the 
search to a single direction, although these techniques 
may be extended to all-sky searches in the future. 

Within this 2500 s x 1000 Hz on-source region, and fol¬ 
lowing (Thrane et al. 2011), we create a spectrogram 
of signal-to-noise ratio SNR(f;/), which is proportional 
to the cross-correlation of the HI and LI strain data. 
Here / refers to the frequency bin in a discrete Fourier 
transform centered on time t. We use ls-long, 50 %- 
overlapping, Hann-windowed data segments, which yield 
spectrograms with a resolution of 0.5 s x 1 Hz. The signal 
to noise ratio is given by 

(SNR(f; /)) = hi (l — cos 2 c) , (56) 

where i is the angle between the direction of propagation 
of the gravitational wave and the plane of the inter¬ 
ferometer (the noise fluctuations can be both positive 
and negative). A formal derivation of SNR(f; /) can be 
found in (Thrane et al. 2011). An example of a SNR(t;f) 
spectrogram with a r-mode injection added to simulated 
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detector noise can be seen in Fig. 10. The problem of 
detection is to identify a track of positive-valued pixels 
in the presence of noise. 

In our simulation, we use a network consisting of two 
4km LIGO observatories assumed to be operating at 
aLIGO design sensitivity; one in Hanford (HI) and one 
in Livingston (LI). We create a spectrogram of signal- 
to-noise ratio SNR(t;f) in our 2500 s x 1000 Hz on-source 
region, which is proportional to the cross-correlation of 
the HI and LI strain data. 
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Fig. 10. — Injection of a r-mode signal using a waveform of f Q = 
1500 Hz and a = 0.1 at a distance of 1 Kpc. The background was 
simulated using MC data colored with the aLIGO sensitivity curve. 


To recover the injected signals, we apply a seedless 
clustering algorithm (Thrane & Coughlin 2013) that 
looks for clusters of positive pixels. The pixels are 
combined to determine the SNR for the entire cluster 
denoted by SNR(c) which is distinct from SNR(t; /) 
associated with individual pixels in a spectrogram. To 
form the SNR for the cluster the pixels are summed 
using their individual sigmas as weights. 

The algorithm was applied to Monte Carlo Gaussian 
noise colored with the aLIGO sensitivity curve, eLIGO 
S5 data (noise) recolored with the aLIGO sensitivity 
curve and Monte Carlo Gaussian noise colored with the 
ET sensitivity curve (Fig.6). The algorithm was used in 
order to determine the threshold for an event with false 
alarm probability < 0.1%. For the MC data colored 
with the aLIGO sensitivity curve we found the cluster 
SNR threshold to be 7.2 while using the ET sensitivity 
curve we found the cluster SNR threshold to be 8.1. 

We then injected simulated signals into the Monte 
Carlo noise at varying distances. We determined the 
distance above which 50% of the signals are recovered, 
which corresponds to the distance at which we can 
observe a signal with false alarm probability < 0.1% 
and false dismissal probability < 50 %. These detection 
distances are shown in table 1 for all 9 different combi¬ 
nations of parameters. 


When testing the sensitivity of a given detection al¬ 
gorithm, we run the algorithm on 10 3 noise maps and 
out of those we set the highest (cluster) SNR as the 
SNR threshold. This threshold corresponds to a fixed 
gravitational-wave strain, h 0 = h t h- Using (14) and this 
threshold hth ., we can test what combination of parame¬ 
ters a and f a and at what distance d the corresponding 
waveform can be detected. Substituting h 0 with hth and 
solving (14) for d we get the distance expression in units 
of Mpc 


where 


d 


1.5 x 1CT 23 


(~ 
\ 1 kHz 



(57) 


/ 3 ~(l + lCT 20 H 2 / o 6 f)-5 (58) 

Equations (57) and (58) suggest that (for a given 
value of h t h) knowing the detection distance d\ for one 
waveform of parameters (aq,/oi) we can estimate the 
expected distance for another waveform of parameters 
(a 2 ,/o 2 ) using 


d 2 


|a 2 | / 1 + AM 2 & y 
M Vl + A|a 2 | 2 / 0 6 J 1 


(59) 


where A ~ 10 20 t. This suggests that when the parame¬ 
ters (a, f 0 ) are such that A|a| 2 / 3 -C 1 then 


d 2 ~ j—4di (60) 

Kl 

This is the case for f a i = 7 x 10 2 , aq = 10 -3 and 
t = 2.5 x 10 3 giving A|ctq| 2 f oX « 2.9 x lO” 6 and 
f o2 = 7 x 10 2 , a 2 = 10" 2 and t = 2.5 x 10 3 giving 
A|oq| 2 /®i « 2.9 x 10 -4 . The seedless based algorithm 
gives for the first case a detection distance of 40 Kpc 
while for the second case it gives 4 Kpc according to (60). 


For A|a| 2 / 3 ~ 1 or A|o| 2 /® > 1 the detection distances 
our algorithm gives follow equation (59). For example 
for f 0 1 = 1.5 x 10 3 , aq = 10 _1 and t = 2.5 x 10 3 giving 
A|oq| 2 /®i « 2.9 and f o2 = 1-5 x 10 3 , a 2 = 10~ 2 and 
t = 2.5 x 10 3 giving A|ck q | 2 /« 3 x 10 -2 we get the 
relation ss 0.19di. The detection distances d\ and 
d 2 are related by d 2 « 0.17di (column 2), d 2 ~ 0.19di 
(column 3) and d 2 ss 0.18di (column 4) of table 1. Thus 
we conclude that the detection algorithm behaves in a 
way consistent with the theoretical predictions. 


12. DISCUSSION 

The most optimistic scenarios in the Bondarescu et al. 
numerical simulations show that r-mode gravitational- 
wave strains take values of 5x 10 -25 Hz -1 / 2 at / = 400 Hz 
and 6 x 10~ 24 Hz -1 / 2 at / = 600 Hz from sources that 
are 100 Kpc away from the detectors. Using (14) we can 
show that these strain amplitudes correspond to r-mode 
saturation amplitudes that range from a = 5.2 x 10 -2 to 
a = 1.9 x 10 -1 . In the worst case scenarios their work 
shows that r-mode gravitational wave strains take values 
of 5 x 10~ 26 Hz -1 / 2 at / = 600 Hz (a = 1.5 x 10" 3 ) 
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TABLE 1 

Sensitivity study using a seedless clustering (SC) algorithm and matched filtering (MF) 


r-modes waveform 
(fo, a) 

(Hz, unitless) 

MC data 
aLIGO s.c./SC 
(Mpc) 

eLIGO data 
aLIGO s.c./SC 
(Mpc) 

MC data 

ET s.c./SC 
(Mpc) 

MC data 
aLIGO s.c./MF 
(Mpc) 

(1500, 0.1) 

1.2 

1.1 

9.7 

22 

(1100, O.l) 

0.97 

1.1 

8.1 

13 

(700, 0.1) 

0.44 

0.42 

4.3 

4.8 

(1500, 0.01) 

0.19 

0.21 

1.8 

2.9 

(1100, 0.01) 

0.13 

0.12 

1.1 

1.3 

(700, 0.01) 

0.040 

0.043 

0.39 

0.54 

(1500, 0.001) 

0.016 

0.021 

0.16 

0.35 

(1100, 0.001) 

0.014 

0.012 

0.11 

0.19 

(700, 0.001) 

0.0040 

0.0044 

0.039 

0.048 


a Distances are calculated for a false alarm probability < 0.1 % and a false dismissal probability <50% using Monte Carlo (MC) aLIGO 
noise and MC ET noise. 

b Doing a background study on 1000 maps the cluster SNR threshold was estimated to be SNR t } l = 7.2 with the aLIGO sensitivity curve 
(s.c.) and SNR t h =8.1 with the ET s.c. Sources are assumed to be optimally oriented. 

c Columns 2, 3 and 4 are the distances at which 50 out of 100 injections were recovered using the SC algorithm and column 5 shows the 
corresponding distances using a highly idealized MF. 


to 2 x 10~ 27 Hz -1 / 2 at / = 350Hz (a = 3.1 x 1CT 4 ). 
Clearly, our results show that we can only be hopeful 
for a large a r-mode gravitational wave detection. 

In the non-linear bulk viscosity mechanism of 
Mark Alford et al. the r-mode amplitudes grow up to 
values of orders 10“ 1 - 1. A 1 KHz gravitational-wave 
signal generated by r-mode oscillations saturated at 
these amplitudes will create gravitational-wave strain 
values from 10 -24 to 10 -23 respectively at distances 
of 1 Mpc where the supernova rate is about 3-4 per 
century. The same wave at a distance of 10 Mpc where 
the supernova rate is at about 1-2 per year (Mannucci 
et al. 2008) will have a strain value of order 10~ 25 to 
10 -24 . Our sensitivity results show that we can exclude 
the 0(1) r-mode oscillations by the absence of any 
r-mode detections at those distances. 


higher than that. 

Initially, this project started with the target to design 
a search that would set an upper bound on the satura¬ 
tion amplitude of the r-mode mass current oscillations 
on neutron stars. For this work we used the Owen et 
al. ’98 model that assumes a polytropic equation of 
state (EOS). Therefore, our work was initially based on 
a specific EOS, using specific values for J 2 and I giving 
Q ~ 0.094. However, we may use different EOS and 
using the integrals for J 2 and I we can calculate their 
corresponding Q value and therefore, the corresponding 
waveforms. Therefore, a hypothetical r-mode gravita¬ 
tional wave detection can not only settle the debate on 
the r-mode oscillations saturation amplitude but also 
impose severe constraints on the EOS of the matter in 
the core of a neutron star. 


Our results show that for a 1 KHz gravitational wave 
generated by a r-mode oscillation saturated at an am¬ 
plitude of order 10 -2 the maximum detection distance 
is 0.1 Mpc. The result for a gravitational wave of the 
same frequency but generated by a r-mode oscillation 
saturated at an amplitude of order 10 _1 is a maximum 
detection distance of 1 Mpc. These results show that 
our search method can detect gravitational-wave strains 
of order down to 10~ 24 . The sensitivity curve of aLIGO 
at 1 KHz is about 5 x 10 -24 Hz -1 / 2 showing that with 
our current algorithms we can detect signals that are 5 
times weaker than the noise. 

Detecting a gravitational-wave strain of order 10” 24 
gives information about the ratio a/d, and does not 
specify the value of a. Therefore, to determine the 
value of a we need an electromagnetic trigger to provide 
the distance for us. Given the distance we can then 
calculate the value of the saturation parameter, a. 
A detection of r-mode gravitational radiation coming 
from a supernova is the best possible outcome of our 
efforts in this search. However, a non-detection can 
also be useful in determining an upper limit on the 
r-mode saturation amplitude and disprove the theories 
that suggest scenarios of r-mode saturation amplitudes 


In the Owen et al. ’98 model, non-linear couplings 
between the r-mode and other inertial modes as well 
as strong magnetic fields (crucial in the evolution/sat¬ 
uration of the r-mode amplitude) were not explicitly 
modeled. However, when compared to the complicated 
numerical simulations in Bondarescu et al. we got a 
significant overlap for the r-mode gravitational wave 
frequency evolution during the first two weeks after the 
neutron star is born. This is the era of a neutron star’s 
life we are interested in. 

An improved model that includes magnetic fields 
effects can be found in (Staff et al. 2012). This work 
considers the same Owen et al. ’98 model but with an 
additional frequency damping term due to magnetic 
braking. The restrictions their findings imply are that 
the r-modes may be dominant only for B-fields of 
magnitude below 10 13 G. This result does not impose 
any changes on the r-mode gravitational waveforms. It 
merely says that we shouldn’t be hopeful to observe 
r-mode gravitational radiation in neutron stars with 
very strong magnetic fields. On the other hand, this 
magnetic braking restriction is dependent on the equa¬ 
tion of state, hence not very decisive on the magnitude 
of the magnetic fields for which we should not hope to 
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detect r-mode gravitational waves. 

The accuracy of the waveforms determines the type 
of decision making algorithm we should be using when 
searching for a signal. Lack of reliable models of 
the physical system implies production of unreliable 
waveforms. This is the reason we did not search for the 
r-mode signal using matched filtering. The algorithms 
we used in this study are based on the statistical 
significance of signal to noise ratios of clusters made of 
pixels above a certain snr threshold. This method did 
not use any knowledge of the signal. Knowledge of the 
r-mode signal can be used and make minor modifications 
in our existing algorithms, however, there was not much 
hope for a dramatic improvement in sensitivity. We were 
still able to recover signals of one order of magnitude 
weaker than the noise. We are currently considering 
algorithms designed to search for specific signals, in our 
case the r-mode. We hope that customizing the search 
algorithm for a specific signal can significantly increase 
the sensitivity in decision making. 

An ideal class of decision making algorithms, which 
are suitable specifically for cases when the signal is 
not precisely (but only crudely) known, is provided by 
machine learning algorithms (MLAs). Three classes of 
MLAs we are currently exploring are artificial neural 
networks (ANN), support vector machines (SVM) 
and local space classifiers (LSC). MLAs are cosidered 
novel methods of signal recognition in the area of 
gravitational-wave searches. We are now working on 
developing r-mode search algorithms using the above 
three techniques. Currently, we are producing data in 
the parameter space of interest. The plan is to repeat 
the same sensitivity study we presented in this paper 
using MLAs for decision making/recognition of the 
desired signals. The results will be compared to our 
present results and published in a future paper. 

13. CONCLUSIONS 

(a) In the worst case scenario of Bondarescu et al. , 
a = 10 -4 — 10~ 3 . If this is correct the results from 
our sensitivity study show that we can only be hopeful 
(with this search method) to detect r-mode gravita¬ 
tional radiation from sources that are up to 21 Kpc away. 

(b) If the best case scenario of Bondarescu et al. (with 
a = 1CT 1 ) prevails then we can be hopeful to detect 
r-mode gravitational radiation from sources that are 
up to 1100 Kpc away. In this case it may not be 
possible to distinguish whether it is Bondarescu et al. or 
Alford et al. mechanism that comes into play. 


(c) To exclude the validity of Bondarescu et al. results 
a r-mode gravitational wave detection yielding a value 
of a > 0.2 is needed. In that case we will have strong 
indications that Alford et al. mechanism is the correct 
r-mode oscillation amplitude saturation mechanism. 

(d) If the best case scenario (of high r-mode saturation 
amplitudes) is correct we may be able to detect r-mode 
gravitational radiation from a supernova event within 
our local group (~ IMpc). Considering the low rate 
of the local group supernovae (3-4 per century) and 
that the latest supernova in the local group occurred in 
January 2014, we may have to wait another 25-30 years 
until the next one. 

(e) The aLIGO and ET sensitivity results show that 
for large amplitudes of a, the detection distances are 
up to 1.1 Mpc (2-3 events per century) and lOMpc (1-2 
events per year) respectively. Using aLIGO, to exclude 
saturation amplitudes of order 10 _1 , a null result at 
a distance of < 0.42 Mpc is required. Using ET, to 
exclude the same saturation amplitudes a null result at 
a distance of < 4.3 Mpc is required. Considering the 
event rates at those distances we can only hope for such 
results only after several years of ET operation. 

(f) We do not need all the complicated physics of 
the neutron star to approximate the r-mode waveforms 
during the first 2 weeks after the onset of the r-mode 
gravitational radiation. A model where high magnetic 
fields are excluded and details of the r-mode saturation 
mechanism are unknown can still give a good approxi¬ 
mation of the r-mode waveforms during the early stages 
of gravitational radiation. 

(g) When the model waveform parameters take values 
over large ranges, choosing a decision making algorithm 
may not be a simple task. Clustering algorithms may 
yield satisfactory results but that is not necessarily 
the best option we can have. Other techniques (for 
example MLAs) that do not require accurate knowledge 
of the waveform parameter values may result in longer 
detection distances. Increasing the detection distance 
by a factor of 2 or 3 can make the difference. 

For lengthy discussions we thank: the Stochastic 
Transient Analysis Multi-detector Pipeline (STAMP) 
group that also developed the code we used, Ruxandra 
Bondarescu, Ira Wasserman, Mark Alford, Lee Lind- 
blorn, Kostas Kokkotas and Nick Stergioulas. 
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